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. Abstract 

We consider an inverse spectral problem for infinite linear mass-spring 
systems with different configurations obtained by changing the first 
mass. We give results on the reconstruction of the system from the 
spectra of two configurations. Necessary and sufficient conditions for 
two real sequences to be the spectra of two modified systems are pro- 
vided. 



Mathematics Subject Classification(2010): 34K29, 47A75, 47B36, 70F17, 
Keywords: Infinite mass-spring system; Jacobi matrices; Two-spectra inverse problem. 
♦Partially supported by CONACYT (Mexico) through grant CB-2008-01-99100 



1. Introduction 



In this work we treat the two spectra inverse problem for Jacobi operators in 
/2(N). The Jacobi operators considered here are obtained from each other by a 
particular kind of rank-two perturbation. The special form of the perturbation has 
a physical motivation; it is the extension to the semi-infinite case of an inverse 
problem for finite mass-spring systems studied in [7j and [18]. 

The Jacobi operator J in the Hilbert space /2(N) is the operator whose matrix 
representation with respect to the canonical basis in Z2(N) is a semi-infinite Jacobi 
matrix of the form 
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where G M and 6„ > for any n G N (see in [2] the definition of the matrix 
representation of an unbounded symmetric operator). J is closed by definition and 
it may be self-adjoint or have deficiency indices (1,1). In this work we deal with 
self-adjoint operators, so, if J 7^ J*, we consider its self-adjoint extensions denoted 
j(9)^ where ^ G M U {00} (see Definition [Da)). If J = J* we assume J^^^ = J for all 
5( G M U {00} (see Definition [T]b)). 

The two spectra inverse problem for Jacobi operators J^^^ takes as input data the 
spectra of two operators in a operator family obtained by perturbing J^^^ in a certain 
way. The solution of the problem is the finding of the matrix (11.11) and the "boundary 
condition at infinity" g if necessary. The case of the operator family consisting of 
rank-one perturbations of a self-adjoint Jacobi operator has been amply studied 
in [51 1121 US] and, in the more general setting of rank-one perturbations of J'-^-' , in 
[ISII23]- Rank-one perturbations can be seen as a change of the "boundary condition 
at the origin" for the corresponding difference equation (see [T9l Appendix]). We 
remark that the case of finite Jacobi matrices has also been thoroughly studied (see 

It is known that the dynamics of a finite mass-spring system is characterized by 
the spectral properties of a finite Jacobi matrix [TT]. Accordingly, in solving the 
inverse problem for mass-spring systems mentioned above, [18] provides necessary 
and sufficient conditions for two point sets to be the spectra of two finite Jacobi 
matrices corresponding to two mass-spring systems, one of which has a mass and 
a spring modified. The results of [I8] are related to the study of microcantilevers 
[2n which are modeled by a spring- mass system whose masses and springs 
constants correspond to the mechanical parameters of the system. The inverse 
problem treated in [18] could be used as a theoretical framework for the problem of 
measuring micromasses with a help of microcantilevers [211 122] . 

Let us consider a semi- infinite spring-mass system with masses {mj}°^]^ and 



1 



spring constants {kj}'jli as in Fig. [H By a standard reasoning (see [HI [IT]) one 
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Figure 1: Semi- infinite mass-spring system 



verifies that the infinite system of Fig. [T] is modeled by the spectral properties of the 
Jacobi operator J with 
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where is the length of the j-th spring at equilibrium. We remark that in [111 IH] 
the obtained matrix corresponds to —J. An alternative physical interpretation is 
provided by a one dimensional harmonic crystal |2H Sec. 1.5]. 

In this work we consider the spectrum of J^^-* to be discrete (if J ^ J* this is 
always the case). Then, the movement of our mechanical system is a superposition 
of harmonic oscillations whose frequencies are the square roots of the modules of 
the eigenvalues. 

Along with the self-adjoint operator J^^^ we consider the family of operators 
J^^\9) {9 > 0) being self-adjoint extensions of the Jacobi operator whose matrix 
representation with respect to the canonical basis in hi^) is 
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J^^\9) (6 > 0) will be the family of perturbed Jacobi operators. Note that the 
operators of the family are not obtained from each other by a rank-one perturbation 
(see ([22D below). 

Going from J^^^ to J^^\9) corresponds to changing the first mass by Am = 
mi{9''^ ~ !)• In other words, 9'^ is the ratio of the original mass mi to the new mass 
mi + Am. This is illustrated in Fig. [2J It is worth mentioning that we also consider 
here the cases when Am < 0, equivalently, ^ > 1, although physical applications 
correspond to ^ < 1 [211 122]. 

The problem of reconstructing the initial and the perturbed matrices by their 
spectra can be then interpreted from the physical point of view as the problem of 
finding the mechanical parameters of the spring-mass system from the frequencies 
of its oscillations before and after the modification. 
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Figure 2: Perturbed semi- infinite mass-spring system 

This work is organized as follows. In Section [2] we lay down the notation, in- 
troduce the Jacobi operators and its perturbations, and present some preparatory 
facts related with the inverse spectral problems of such operators. Section [3] gives 
an account of the spectral properties of the family of perturbed Jacobi operators 
J^^\6). The problem of reconstruction is treated in Section |H This section gives 
some necessary conditions for the spectra of J^^\9), provides an algorithm for re- 
construction of the matrix and establishes uniqueness of the reconstruction. Finally, 
Section [5] gives necessary and sufficient conditions for two sequences of real numbers 
to be the spectra of J^^^ and its perturbation j'^^\0) {6 ^ 1). 

2. Preliminaries 



Let T be a second order symmetric difference expression such that for any se- 
quence f = {fk]ken 



(T/)i 



qifi + , 

bk-lfk-l + Qkfk + bkfk+1 



keN\{l}, 



(2.1) 
(2.2) 



where, for n G N, 6„ is positive and g„ is real. Let /fin(N) be the linear space of 
complex sequences with a finite number of non-zero elements. In the Hilbert space 
/2(N), let us consider the operator whose domain is /fin(N) and acts as the expression 
T. This operator is symmetric since it is densely defined and Hermitian, and thus 
it is closable. Now, let J be the closure of this operator. 

We have defined the operator J so that the semi-infinite Jacobi matrix f 1 1.1 1) is 
its matrix representation with respect to the canonical basis {6n}'^=i in hi^) (see ^ 
Sec. 47] for the definition of the matrix representation of an unbounded symmetric 
operator). Indeed, J is the minimal closed symmetric operator satisfying 



{JSn, Sn+k) = {^n, J ^n+k) = , 



nGM, A;GN\{1}. 



We shall refer to J as the Jacobi operator and to (11. ip as its associated matrix. 
The operator J* turns out to be given by 



dom( J*) = {/ G : T/ G UN)} rf = Tf 
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which follows directly from the definition of J [T, Chap. 4 Sec. 1.1], [201 Thm. 2.7]. 
If one gives the complex number /i, the solution of the difference equation, 

(T/) = c/, CgC, 

is uniquely determined from (12. ip and (12. 2 p by recurrence. For the elements of this 
solution when /i = 1, the following notation is standard [11 Chap. 1, Sec. 2.1] 

P„_i(C):=/„, neN, 

where the polynomial Pk{C) (of degree k) is referred to as the k-th orthogonal poly- 
nomial of the first kind associated with the matrix (II. ip . The polynomials of the 
second kind QkiC) are defined as the solutions of 

(T/)fc = C/fc fcGN\{l} 

under the assumption that /i = and /2 = b^^. 

The sequence P{() := {Pk-i{C)}T=i ^finll^), but it may happen that 

oo 

5^|Pfc(C)|'<oo, (2.3) 

fc=0 

in which case P{C) G ker(J* — (I). Since J is symmetric, if the series in (12.30 is 
convergent for one ( in the upper half plane C+ (the lower half plane C_ ) , then it is 
convergent in all C+ (C_). Actually, because of the reality of Pfc-i(C) ^^i all k eN, 
the series in (12.31) will then be convergent in all C \ M and J has deficiency indices 
(1, 1). When the series in (12.31) is divergent for one C in C\M, J has deficiency indices 
(0, 0) and the operator is self-adjoint since J is closed. There are known conditions 
on the matrix (II. ip which guarantee that J is self-adjoint [U Chap. 1, Addenda and 
Problems], [3', Chap. 7, Thms. 1.2-1.4]. 

We now introduce the operators that will be at the center of our considerations 
in this work. 

Definition 1. Let the operator J^^^ be defined as follows: 
a) In case J 7^ J*, define the sequence v{g) = {vk{g)}kLi such that Wk eN 

vk{g) := Pfc-i(O) + gQk~i{0) , geR 

and 

Vk{oo) := Qk-i{0) . 
Let J^^^ be the restriction of J* to the set 

i / = {fk}ken G dom(J*) : lim h{vk{g)fk+i - fkVk+i{g)) = [ . 
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When ^ G M U {oo}, J'^s) 

runs over all self-adjoint extensions of J. Moreover, 
different values of g imply different self-adjoint extensions [211 Lemma 2.20]. 

b) In case J = J*, let us define J^^^ := J for all G M U {oo}. 

Alongside the operator J^^^ we consider the operators Jn^ {n G N) in the Hilbert 
space /2(N) 0span{(5i, . . . , 5„} defined by restricting J^^'> to hi^) 0span{5i, . . . , 5^}. 
Thus, J!f^ is a self-adjoint extension of the Jacobi operator whose associated matrix 
is (11.11) with the first n columns and n rows removed. 

Finally we introduce the perturbed operators J*^^^ (9). They are defined as follows. 
Consider J^^^ with fixed G M U {c>o} and take any ^ > 0. Then 

J(a)^e) := J(3) + qiiO' - 1) ■) 6, + hie- l)((5i, ■) 62 + {62, ■) 5i) , (2.4) 

where we take the inner product to be antilinear in its first argument. By this 
definition j(9)(0) is a self-adjoint extension of the Jacobi operator whose associated 
matrix is ([L2l). Note that j'-s^O) is a finite-rank perturbation of J'-^^ and thus 
dom(j(9)) = dom(j(^')(^)). 

Fix G M U {00} and take the resolution of the identity E^^\t) of J^^\ so 

= [ tdE^3\t). 



Since J^^^ is simple [H Sec. 2.2 Chap. 4], it is particularly useful to consider the 
function 

p(3)(t) := (5i,£;(5)(t)5i) , tGR. (2.5) 

It turns out that all the moments of the measure generated by p^^^ are finite [H 
Thm. 4.1.3], that is, 

Sk= f t^dp^^\t) < 00 VA; G N U {0} , (2.6) 

and the polynomials are dense in L2(M,rfp(3)) lli Thms. 2.3.2, 4.1.4], [20, Prop. 4.15]. 
In this work we also make use of the so-called Weyl m-function 



m 



{9) 



(C):=(5i,(J(^)-C/)-'5i> , C^^(>^). (2.7) 



The functions (12. 5p and (12.71) are related by the Borel transform, viz.. 



m(^)(C) 

so rrS^'> is a Herglotz function, i.e.. 



dfMHf)_ 



> , Im C > 

ImC 



Using the von Neumann expansion for the resolvent (cf.pH Chap. 6, Sec. 6.1]) 

k=o ^ ^ 

where C G C \ a( one can easily obtain the following asymptotic formula 

m^H0 = -l-f,-^^ + OiC'), (2.8) 

as C — ^ oo (Im^ > e, e > 0). 

The inverse Stieltjes transform allows to recover the spectral function fl2.5p from 
its corresponding Weyl m-function (12. 7p . So they are in one-to-one correspondence. 
Furthermore, either (12. 5p or (12. 7p uniquely determines the Jacobi operator J^^\ i.e., 
the matrix (II. ip and the parameter g in the non-self- adjoint case. Indeed, there 
are two general methods for recovering the matrix (II. ip that work without any 
assumption on the spectrum. One method, developed in [3] (see also |23]), makes 
use of the asymptotic behavior of the Weyl m-function and the Riccati equation [HI 
Eq.2.15], [231 Eq.2.23], 

6>(f)(C) = gn-C-^J— -, neN, (2.9) 

where mn\Q is the Weyl m-function of the Jacobi operator Jn^ (mg = m). 

The other method of reconstruction (see [3l Chap. 7, Sec. 1.5, particularly, Thm. 1.11]) 
has its starting point in the sequence {t'^j^Q, t G M. From (12. 6p all the ele- 
ments of the sequence {t'^j^Q are in L2(R,dp^9)^ and one can apply, in this Hilbert 
space, the Gram-Schmidt procedure of orthonormalization to the sequence {t'^j^Q. 
One, thus, obtains a sequence of polynomials {Pfc(^)}fclo normalized and orthogo- 
nal in L2{M., dp^^^). These polynomials satisfy a three term recurrence equation [3l 
Chap. 7, Sec. 1.5], [201 Sec. 1] 

tPk-i{t) = bk-iPk-2{t) + qkPk-i{t) + hPk{t) , A;GN\{1}, (2.10) 
tPo{t) = qiPoit) + hPiit) , (2.11) 

where all the coefficients bk {k G N) turn out to be positive and qu {k G N) are real 
numbers. The system (I2.10p and (12. lip defines a Jacobi matrix which is either the 
matrix representation of either J^^^ or a restriction of J*^^^ depending on whether 

J ^r. 

The function (12. 7p . equivalently (12.51) . determines the parameter g which defines 
the self-adjoint extension when the reconstructed matrix turns out to be the matrix 
representation of a non-self-adjoint operator. Indeed, consider a pole 7 of rrS^^ (there 
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is always one when J ^ J*) and evaluate Pk{l), A; G N. Then either 



lim hk{Qk-MPk{l) - Pk^i{l)Qk^M) = 

k—^OD 

which means that (7 = oo, or 

\imk^^bkiPk-iiO)Pkh) - Pk-i{l)Pk-im 
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limfc^oo bkiQk-MPkii) - Pk-ih)Qk-iiO)) 



The details of this recipe are explained for instance in Sec. 2]. 

Since any simple self-adjoint operator in a infinite dimensional Hilbert space is 
unitarily equivalent to some operator J = J* [1| Thm. 4. 2. 3], [2] Sec. 69], in the case 
J = J*, aiJ'^a)) 

may be any non-empty closed infinite set in R. In particular J^^^ 
may have discrete spectrum, that is, cress (J*-^-*) = 0. When J ^ J* ^ this is always 
the case, that is all self-adjoint extensions J^^^ of the non-self-adjoint operator J 
have discrete spectrum [231 Lem. 2.19]. 

Assume that J has discrete spectrum (this always happen if J 7^ J*), so the 
spectrum is a sequence of real numbers, {Afc}fc, without finite points of accumulation. 
The simplicity of J'-^-' implies that all eigenvalues are of multiplicity one. In this 
case the function p^^\t), defined by (12.51) . can be written as follows 

P^'\t) = J2-^ (2.12) 

where the coefficients {ak}k are called the normalizing constants and according to 
[U Chap. 7, Thm. 1.17] are given by 



00 

k=0 



^ = J2\Pk{Xn)f . (2.13) 



Thus, from f l212|) and one has that 

'(C) = E (2.14) 
^ ak{\k - C) 



Remark 1. In the case of discrete spectrum, the set of poles of the meromorphic 
Weyl m-function coincides with a{J^9^). By (12. 9p . the set of zeros coincides with 
cr(j|^^). The zeros and poles of the Weyl m-function are simple and interlace as 
occurred to any meromorphic Herglotz function. Interlacing means that between 
two contiguous poles there is exactly one zero and between two contiguous zeros 
there is exactly one pole (see the proof of ^LQ\ Chap. 7, Sec. 1, Thm. 1]. 

Remark 2. By elementary perturbation theory (Weyl theorem), J^^^ has discrete 
spectrum if and only if J^^^ (6) has discrete spectrum. Note that J^^^ (6) has simple 
spectrum since it is a self-adjoint extension of a Jacobi operator. 
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3. Direct spectral analysis of J^^^ and J^^\9) 



We begin this section by noting that 

= j[9){^d)^ V^>0. 

Fix 5f G M U {00} and consider the Weyl m-functions m^^\ m^^^ of the operators 
J*^^^ and J^^\6). Therefore, taking into account that mf'' and mf'' coincide, (12 .Qp 
imphes that 

1 



Let us now consider the function 

Remark 3. In view of Remark [21 if J*^^^ has discrete spectrum, the function m is 
meromorphic by (13.21) . Since the zeros of m*^^) and m*^^) are the same (see Remark[T]), 
it follows that for all 6^ > the set of poles of m is a subset of <j( J*^^)), while cr( J*^^^ {9)) 
contains all the zeros of m. Observe also that, from (13.11) . G a{J^^^) if and only if 
G a{J^^\6)). Moreover, whenever ^ 7^ 1, (13. ip implies that the sets a{J^^^) and 
a{J^^\6)) can intersect only at 0. 

Remark 4. By [151 Chap. 7, Sec. 3.5, Thm. 3.9] the zeros of m are analytic functions 
of the parameter 9. The same is true for the eigenvectors of J^^\6). 

Proposition 3.1. Let J^^^ have discrete spectrum and let {XkiO)}k be the set of 
eigenvalues of J^^\6) (9 > 0). For a fixed k the following holds 

where ak{0) is the normalizing constant corresponding to Xk{6)- 

Proof. Let us denote by f{9) the eigenvector of J^^\6) corresponding to Xk{0)- We 
assume that f{d) is normalized in such a way that 

(5i,/W) = l. (3.3) 

Pick any small real r (it suffices that |r| < 9). Then, taking into account that 
dom(J'^^)) = dom( J*^^)(^)) and the self-adjointness of J^^\6) for any ^ > 0, we have 
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that 

{X,{9 + r) - A,(^)) {f{e), f{9 + r)) = {f{e), J^^\e + r)f{6 + r)> 

-{j^^){e)f{e)j{e + T)) 

= (fie), ij^a\e + r) - + J^^\e))fi9 + t)> 

= {f(9),iJ^^\9 + T)-J^^\9))fi9 + T)) . 

From (13.31) it follows that the entries f{9 + T) and f{9) are the polynomials of the 
first kind associated to the matrix of J^^\9 + r) and J^^\9), so 



f2{9 + T) 



Xki9 + T) 



f2{e) 



\k{9) - 9\, 



{9 + T)h " ' 9bi 

Now, taking into account these last equalities and (13. 3p . together with 



f{29T + T^)qi rbi • ■ A 



J^3\9 + t) - J^'^\9) 














V 
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one obtains that 

{X,{9 + r)-X,{9)){f{9),f{9 + T))- 

Therefore, on the basis of Remark HI one has 

X,i9 + r)-X,i9) ^ 1 

r^o{f(9)Ji9 + r)) 



Xk{9 + T) ^ Xk{9) 



Xk{9 + r) , Xk{9)\ 2Xk{9) 



9 + T 
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9aki9) 



□ 



The proposition below can be proven by means of Remark [3], HI and Proposi- 
tion [3lT] However, we present an alternative proof based on the following expression 

miO = - l W'\0 + e\ (3.4) 
which follows from (13. ip and (13. 2p . 

Proposition 3.2. Fix g E MUjcxo} and let J^^^ have discrete spectrum. The spectra 
(r{J^3)^, a{J^^\9)) interlace in M+ and M_. Moreover, a(j(s)(^)) in M+ (M_; is 
shifted with respect to a{J^^^) to the left (right) if 9 < 1, and to the right (left) if 
9>1. 
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Proof. In view of Remark [3l one only needs to verify that between two positive 
and contiguous eigenvalues of J^^^ there is only one eigenvalue of J^^\6). Take two 
positive and contiguous eig envalues of A < A. Due to fimj) . one has 

liin m^^\t) = +00 , lim m'^^^t) = -00 . (3.5) 

Now, in fl3.4p assume that 6 > 1. Thus, because of the positivity of A, A, fl3.4p and 
(13.51) imply that 

lim m{t) = +00 , lim m{t) = —00 . 

tm teR 

Since m is analytic on the interval (A, A), it should cross the 0-axis an odd number of 
times. If it crosses this axis three or more times as in Fig. [3] (a), then, by Remarks [1] 
and[3], there are at least two elements of a{j[^^) in (A, A). But, because of RemarklH 
this would contradict the fact that A, A are contiguous. Observe that one should 






Figure 3: Impossible crossings of the 0-axis by m 

discard the possibility of one crossing of the 0-axis and a tangential touch of it as in 
Fig. [3] (6) and (c). But again the impossibility of this follows from the fact that the 
poles of m^^^ are simple (see Remark [1]). By the same token the spectra interlace 
in M_. The case ^ < 1 is treated in a similar way. The second assertion follows 
directly from Proposition 13.11 □ 

Remark 5. We note that a{J^^^) fl R+, a{J^^^) fl M_, may be finite or empty. 

4. Inverse spectral analysis for J^^^ and J^^\6) 

In this section we find some necessary conditions for the spectra of J^^\6) [6 > 
0). Also we provide a reconstruction algorithm of the Jacobi matrix and establish 
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uniqueness of the reconstruction. Some of the formulae obtained in this section (see 
for instance Corollary 14. ip have an analogous one in the finite case [7], [18] . 

A central part of our approach is the Weyl m-function and its properties. We 
begin our discussion by setting out a convention for enumerating the elements of 
the spectra. 

Convention. For a given countable set of real numbers S without finite points of 
accumulation, let M be an infinite subset of consecutive integers such that there is 
a one-to-one increasing function h : M ^ S with the property that, h~^{0) = 0. 
Thus, M is semi-bounded from above (below) if and only if the same holds for S. 
We write S = {Xk}keM, where = h{k). Note that in the sequence {Xk}keM only 
Ao is allowed to be zero. Thus, if —1,1 G M, then 

A_i < < Ai . 

In the sequel, the spectra of all operators will be enumerated according to this 
convention. 

When {Xk}keM is considered together with a sequence interlacing with it, we 
use the same set M for enumerating both sequences. For instance, if {Xk}keM and 
{lJ,k}kGM are interlacing and not semi-bounded, then one can assume that 

Xk < fik < Xk+i , Wk e M. 

The following auxiliary result can be found in Sec. 4]. We sketch the proof 
here for the reader's convenience. 

Lemma 4.1. Let J^^^ have discrete spectrum and assume that a{J^^^) = {Xk}keM, 
and a{j[^^) = {rik}k&M- Then, the following formula holds for the Weyl m-function 
of J(3) 



-1 



Moreover, C < and 

Vk<Xk< r]k+i WkeM, (4.2) 
if a(J^^^) is semi-bounded from above, while, C > and 

Xk<Vk< Xk+i WkeM (4.3) 

otherwise. 

Proof. Assume first that a{J'^^^) is semi-bounded from below. Since the greatest 
lower bound of J does not exceed the greatest lower bound of the smallest 
element of {Xk}keM is less than the smallest of {rik}keM (see [4, Chap. 6 Sec. 1.3]). 
Thus one can enumerate the sequences {Xk}keM and {r]k}k&M so that they obey our 
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convention and (14. 3p . According to [16, Chap. 7, Sec. 1, Thm. 1], (14.11) holds with 
C > 0. 

Clearly, when a{J^^^) is not semi-bounded, the sequences can be arranged to 
obey (53]), and then (gl]) holds with C > 0. 

Now suppose that is semi-bounded from above. Then a{—J^^^) is semi- 

bounded from below and, consequently, the greatest of {rik}keM is less than the 
greatest of {\k\keM- Thus {AfcjfcgM, and {r^fcjfceAf cannot be arranged according 
to fl4.3p . However, we are still able to use fl4.3p for arranging the zeros and poles 
of the meromorphic Herglotz function — that is, we use (14. 2p . Therefore [T6l 
Chap. 7, Sec.l, Thm. 1] gives 



k^O 

For completing the proof it only remains to note that the last equation can be 
rewritten as asserted in the lemma. The infinite product in (14. ip is convergent 
because of (14.20 (see the proof of [16, Chap. 7, Sec.l, Thm. 1]). □ 

Another auxiliary simple result to be use later is the following lemma. 

Lemma 4.2. Let J^^^ have discrete spectrum and {\k{Q)}k be the set of eigenvalues 
of J^s\d). Then, the series 

Y ^ (4.4) 

keM '^^ ' 

converges uniformly in [^i,^^2] C M+ to Si{d) (see Ii2.6\) ). 

Proof. From (12. 6p and (12.121) . it follows that the series converges pointwise to si{9). 
The series 

k&M ' 

converges also pointwise to the function s^iff)- Since this function is continuous in 
[6'i, 6^2], then (14. 5 p is uniformly convergent in that interval (see ^25, Sec. 1.31]). Now, 
for any Q G [6'i,6'2] and |Aa;| > 1, one has 

|Afc| < , 

so (14.40 is uniformly convergent in [^?i,02]- D 

Remark 6. Proposition 13.21 tells that the interlacing of the sequences a{J^^^) = 
{Xk}k and a{J{9)) = {fik}k is different in M+ and M-. So let us agree to enumerate 
the sequences according to our convention (the subscripts of the sequences run over 
M and only the eigenvalues with subscript equal zero are allowed to be zero) and 
obeying 

Xk < ^k < Afe+i in M+ , fXk < Xk < fJ-k+i in M_ , 
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when 6 > 1, and 

fJ'k < < /Jfc+i in IR+ , < l^k < Afc+i in M_ , 

if 6* < 1. 

Proposition 4.1. Fix g E W U {oo} and < < 62- Let J^^^ have discrete 
spectrum and assume that a{J^^\6i)) = {Xk}k(^M and a{J^^\92)) = {/UfcjfceM, where 
the sequences have been arranged according to Remark{^ Then, 



J2if^k - Xk) = qi{0l - Ol) . 



keM 

Proof. Observe that from Proposition 13.11 it follows that 



fJ-k — Xk 



ei 



Consider a sequence {M„}5^]^ of subsets of M, such that M„ C Mn+i and U„M„ 
M. Thus 

Xk{0) \ d9 



5^(/ifc-Afe) = 2 hm [\Y1 

keM ''^^ \fceM„ 



ak{e) 9 



By Lemma 14.21 and the fact that 

s,{9) = {d^,J^'\9)6,)=q^9^ 

one obtains 

f02 



- Xk) = 2qi r 9d9 = qi{9l - 9l) 
keM "^^1 



□ 

rum 



Proposition 4.2. Fix g eMU {00} and < 9 ^ 1. Let J'^^^ have discrete spect 
and assume that a{J^^^) = {Xk}keM ^Li^d a{J^^^{9)) = {fik}keM, where the sequences 
have been arranged according to Remark\^ Then, 

keM ^ 

Proof. Consider a sequence {M„}J^^ of subsets of M, such that M„ C Mn+i and 
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U„M„ = M. From (jH]) and ([321) it follows that 



m(C) = C- — lim 



'"•■""n(-a('-a" 

On the other hand, by Proposition I4.H it holds true that 

k^O kj^O k^O 

From (El]) and it follows that 

lim m(C) = l. (4.8) 

ImC>e>0 

Also, on the basis that the second product on the r. h. s of (14 .7^ converges uniformly, 
one has 

lim TT^= lim TT + =1- (4-9) 

ImC>efceA^ ImC>£'=eM \ 

Thus, dH]), (i2D, dlS]), and (gS]) imply that 

fcGM 

and the proposition is proven. □ 

Corollary 4.1. Fix g E MU {oo} and 9 > 0. Let J^^^ have discrete spectrum and 
assume that a{J^^^) = {Xk}k o.nd a{J^^\9)) = {iJ.k}k, where the sequences have been 
arranged according to Remark\^ Then, 



V- fJ-k 



ri — \u 

keM ' 
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where rj is any element of a{j[^^). Moreover, when ^ a{J^^^), 



e^=\{^- (4.10) 



fceM 



and, ifOe ct(J(»)), 

^^ = -^<'ao II , (4.11) 



keM 



where is given in Ii2.13\) . 



Proof. The first two identities for 6^ are a straiglitforward consequence of Proposi- 
tion |1]2] and (13. 4p . As regards to fl4.1ip . note that, from f l2.14p . one has 



a,^ = -Resm(C). (4.12) 



Thus, according to (13. 4p . 



^^-«o^(^^-l) = m(0)= n Y^- (4.13) 

□ 

Remark 7. Due to (I4.13P and the properties of the normahzing constants, when 
one the following inequalities hold depending on the value of ^ 7^ 1: 



9^ < m(0) = 17 ^ < 1 1< m(0) = T] ^ < 9^ 

keM " k£M 

k^O k^O 



Theorem 4.1. Fix g E M.U {00} and 6 > 0. Let J^^^ have discrete spectrum and 
assume that0^a{J). The spectra a {J^3)^, ct(>)(^)) (9^ 1) uniquely determine the 
Jacobi matrix ( (i.ip . that is the operator J, the parameter 9 defining the perturbation, 
and the parameter g specifying the self-adjoint extension when J ^ J* . 

Proof. Given the sequences a{J^^^) and a{J^^\9)), one finds the parameter 9 from 
(I4.10p . Proposition 14.21 yields the function m and equation (13. 4p the Weyl function 
m^^\ According to the Preliminaries this function allows to recover the matrix 
associated to the Jacobi operator and the parameter g which determines the self- 
adjoint extension when J ^ J*. □ 

Theorem 4.2. Fix g E {00} and 9 > 0. Let J^^^ have discrete spectrum 
and assume that E a{J). The spectra a{J^3)^, a{J^<^\9)) (9^1), together with 
either qi or a^, uniquely determine the matrix associated to J, the parameter 9, 
and the parameter g when J ^ J* . Alternatively, the spectra a{J^^^), a{J^^\9)) 
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and the parameter 6^1 uniquely determine the matrix corresponding to J and the 
parameter g if necessary. 

Proof. This follows immediately from the proof of the previous theorem, taking into 
account (14. lip . Note that 9 can be determined either by Proposition 14. II or by the 
asymptotic formula 

m(C) = l + «i<l^ + 0(C-), 
as C — ^ cxD (ImC > e, e > 0), obtained by combining (12 .Sp and (13 ■4p . □ 



5. Necessary and sufficient conditions for the spectra 

of J^s) and 

The following statement gives an if-and-only-if criterion for two sequences to be 
the spectra of J^^^ and j'^^\9). In the finite case the interlacing condition given in 
a) (see below) is necessary and sufficient [7j,[T8]. 



Theorem 5.1. Given two infinite real sequences {\k}k o-nd {/ifcjfc without finite 
points of accumulation, such that none of them contains the zero, there is a unique 
positive 6, a unique operator J, and a unique g G MUjoo} if J ^ J* , such that 
{/^fc}fc is the spectrum of J^^\9) and {Xk}k is the spectrum of J^^^ if and only if the 
following conditions are satisfied. 

a) {Xk}k and {fj.k}k interlace in M+, M_ with one sequence shifted to the right (left) 
in ]R+, (M_j with respect to the other one. Thus, the sequences can be ordered 
according to Remark\^ 

b) The following series converges 

fcsM 



By condition^) the products TT y~ T~"; TT T~ are convergent, so define 

At — An Ah 

k&M k€M 
k^n 



ifJ'n - A„) Y\ 



keM " 

^^"^ -, yneM. (5.1) 



\keM 
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c) The sequence {r^jngM is such that, for m = 0, 1, 2, . . . , the series 

Tfc converges. 

keM 



d) If a sequence of complex numbers {/3k}k&M is such that the series 



kf Tk converges 



kGM 

and, for m = 0,1,2, ... , 

J2 ^'^>^kTk = , 

fceM 

then = for all keM. 

Proof. In view of Propositions 13 . 2 1 and 14 . ![ for proving the necessity of the conditions, 
it only remains to show that for all n G M, Tn = Indeed [cj) and[^ will follow 
from the fact that all moments of the spectral measure (12.121) exist and that the 
polynomials are dense in L2(M, p). 

From (Q, (E2]), and Proposition O , it follows that 



n 



An (6*^ - 1) ^J^j An - Afc 

Hence, by Corollary 14.11 one verifies that = a~^. 

We now prove that conditions [oj) , E]) , [cj) , and[^ are sufficient. 
The condition [aj) implies that 

^^^^>0 ykeM,k^n 

An ~ Afc 

On the other hand, by \E) one can define the number 



k€M 



which is clearly positive and also k > 1 ii \^k\ > \W for all A; G M and k < 1 ii 
< |Afc| for all k G M. Thus, 

>0 VnGM 



\Jk - I) 
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Hence, for all n G M, Tn > 0, so define the function 

Pit) ■.= J2n. (5.3) 



Xk<t 

It follows from [cj) that the moments of the measure corresponding to p are finite. 
Now, on the basis ofdj) and El), define the meromorphic functions 

s(c) n ^ 

and 

ffi(c) - n £ 

ff,(^Q,= . . (5.4) 



Thus, taking into account (15. ip . one has 

Res m(C) = f n f - ^^t^T(C) = -r„ . (5.5) 



\keM 

In view of what was done earlier 



lim m(C) = 1 . (5.6) 

(— s>oo 
ImC>e>0 



Therefore, 



lim m{0=lY\^-l] lim ^ = (5.7) 

ImC>e>0 \keM / ImC>e>0 



By (15. 5p and ( 15. 7p . [16, Chap. VII, Sec.l Theorem 2] implies that 

n 

On the other hand, using (15.61) . one obtains 

ImC>e>0 " / ImC>e>0 V ^SM 

But 

lim Cm(C) = - ^fc , 

ImC>e>0 'seA^ 



so it has been proven that, for the function given in fl5.3p . 



/ dpit) = 1 . 
Jr 



Thus the measure corresponding to p is appropriately normahzed and all the mo- 
ments exist, so in L2(M, p) apply the Gram-Schmidt procedure of orthonormalization 
to the sequence {t^j^Q to obtain a Jacobi matrix as was explained in the Prelimi- 
naries. Denote by J the operator whose matrix representation is the obtained matrix 
(cf. pi Sec. 47]). Now, depending on the sequence of moments, J is self-adjoint or 
not. If J = J*, the function p is the resolution of the identity of J, while if J 7^ J*, 
p corresponds to the resolution of the identity of a self-adjoint extension of J. This 
is a consequence of condition since it means that the polynomials are dense in 
L2(M,p) [i Prop. 4.15]. 

Finally, denote by J^^^ the self-adjoint extension of J corresponding to p and 
consider the operator j'^^\6) obtained from J*^^) as indicated in the Preliminaries 
with 9 given by fl4.10p . By construction the sequence {Xk}k(^M is the spectrum of 
Por the proof to be complete it only remains to show that {pk}k&M is the 
spectrum of J^^\9). For the function given in f l3.2p . taking into account (13.41) and 
(^JS>, one has 

0' + C{O'-l)Yl ^ 



m(C) 



akih - C) 



On the other hand, from (15. 4 p and (15. 8p . it follows that 



keM 



Xk-C 



But we have already proven that ^ = for k G M. Thus m 
the zeros of m are given by the sequence {pk}keM- 



m, meaning that 
□ 



Theorem 5.2. Let {Xk}k o,nd {pk}k be two infinite real sequences without finite 
points of accumulation, such that each of them contains exactly one element equal 
zero, and consider any positive real number 6^1. There exists a unique operator 
J, and a unique (7 G M U {00} if J ^ J* , such that {pk}k is the spectrum of J^^\9) 
and {Xk}k is the spectrum of J^^'' if and only if the conditions \^) , \^) , and\^ hold 
with 

Pn — Xn YJ Pk ~ Xn . 

- ■ 11"; ' neM,nf^O. 



A„(6'2-l) 11 Xk-X 

" ^ ^ keM " 

ky^n 



Xk 



keM 
k=iO 
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where 



^ TT T~ ^-^ {/^fc}fc ^5 shifted to the left 



> 



fcgM 

n 

fcGM 



w.r.t. {Xk}k 



otherwise. 



(5.9) 



Proof. The proof is analogous to the proof of Theorem 15. 1[ Recall that by our 
convention for enumerating the sequences Xq = fiQ = 0. Thus, for proving the 

^ in^ and (Eal) 



necessity of the conditions [ij)-^), one only should verify that tq 
holds. This is immediate in view of fl4.13p and Remark [71 The sufficiency of the 
conditions is established as in the proof of Theorem 15. II Here, one substitutes f l5.2p 
by 

■Q f^k 



K 



keM 



A. 



and (15^ by 



m(C) 



m(C) - 0' 



Then, one verifies that ReS|j=A„ ^(C) = ~'^n for all n G M and J2keM'''k = 1- Note 
that ( 15. 9p guarantees that r„ > for all n G M. The rest of the proof repeats 
that of Theorem 15.11 taking into account that now the zeros of m are given by 
MkeM \ {0}. □ 

Theorem 5.3. Given two infinite real sequences {Xk}k CLnd {^k}k without finite 
points of accumulation, such that none of them contains the zero, there is a unique 
positive 6 and a unique operator J = J* such that {jj.k}k is the spectrum of J^^\9) 
and {Xk}k is the spectrum of J if and only if conditions [^), [^), together with 



det 



/ So 

Sl 



Sl 
S2 



Sn \ 

Sn+1 



d') 



lim 

n—>oo 



\Sn Sn+1 



S2n J 



det 



/ S4 

S5 



55 
56 



Sn+2\ 

Sn+3 



0, 



\5n+2 5„4.3 ■ ■ ■ S2n / 

where s„ := ^^.^^ X^t^ for n in NU {0} are fulfilled. Note that by our convention 
on the notation J^^^O) is a non-singular finite-rank perturbation of J which does 
not depend on g. 

Proof. We again repeat the reasoning of the proof of Theorem 15. 1[ Clearly, Sn {n G 
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NU {0}) are the numbers given in (12.61) . Thus, on the basis of Hamburger criterion 
(see [U Chap. 2 Addenda and Problems, Example 9]), d') holds when J = J*. For 
the sufficiency, note that, due to Chap. 2 Addenda and Problems, Example 9], d') 
implies that the measure corresponding to the function given in (15. 3 p is the unique 
solution of the moment problem, so J = J* and d) is not needed. □ 

Remark 8. Admittedly, d') is not easy to check, however it allows to give necessary 
and sufficient conditions in the self-adjoint case. Note that one can also give the 
analogous self-adjoint version of Theorem 15.21 by substituting condition d) for d'). 
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